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(1) Introduction 



This software package generates the element connectivity and nodal 
point coordinates for two dimensional isoparametric 8 nodal point 
finite elements. 

The maximum number of finite elements is limited only by the 
available disk space. 

Furthermore several techniques such as 

- nested menus 

- keyboard filtering 

- graphics support 

- preparation of screen output in RAM memory 

- dynamic memory allocation 

- module communication via disk files 

are employed in order to achieve a user friendly, fast executing 
system that takes full advantage of small machines with limited 
RAM memory. 



(2) Hardware requirements 

The hardware configuration necessary to run this program is an 
IBM PC or compatible with 256K of RAM memory and 2 DS/DD disk 
drives. The PC must be equipped with a video card that supports ! 
the high resolution graphics mode. Advantageous but not necessary | 
is a hard disk. ' 

Depending upon what floating point library is linked in with the i 
compiled code, the program may be used with the 3087 (rsp. 80287 
for the PC/AT) math coprocessor chip or an 8086/8088 (rsp. 80286 
for the PC/AT) chip. 



(3) The mesh generator concept 

The theory of isoparametric elements is discussed in full detail 
in ref.Cl] (pp.9-76). We will use the notation of ref.[l]. 

An example shall demons t ratethemesh generator concept. 

Let us assume we want to discretize the region shown in Fig. 3.1. 

Fig .3.1; y 



I 
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(3.1) Step 1 : 



"Split" R into cartesian super elements (CSEL). 



Each CSEL corresponds to one isoparametric quadratic 3 nodal point 
real element (see ref. C 1 ]> P • 21/39/104 ) . 

We could choose two CSELs as indicated in Fig. 3.2. 



Fig. 3.2 




B 



A, PI , P2 , P3 , B are 
located on the 
circle boundary 



The circle boundary AB in Fig. 3.1 is now approximated by the two 
quadratic polynomials passing through (A,P1,P2) and (P2,P3,B). 

We get a "better" approximation if we choose 4 CSELs . 



Fig . 3.3 




B 



A, PI , . . , P7 , B are 
located on the 
circle boundary 



Now the circle boundary AD is approximated by the 4 quadratic 
polynomials defined by ( A , P 1 , P 2 ) , ( P 2 , P 3 , P4 ) , ( P 4 , P 5 , P 6 ) and 
(P6,P7,B) . 

It should be mentioned that in both cases (2 CSEL/Fig.3.2 and 
4 CSEL/Fig. 3 . 3 ) the remaining boundary BCDE of R is represented 
exactly by the Snodal point elements. 

Depending on how "accurate" we want to approximate the actual 
boundary of the given region R, we have to "split" R in adequately 
many CSEL ’ s . 

If R consists of different materials, we have to decompose R3 with 
respect to the material and "split" each piece (now consisting of 
one material) into CSELs. 
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(3.2) Step 2 : 



Transform the CSELs into master grid super 
elements (MSEL). 



Themaster grid consists - in this program - of 10x10 MSELs. Each 
MSEL has a unique identification number as shown in Fig. 3.4 : 



Fig. 3.4 ; 



identification numbers of the MSELs 



global 

c 

0 

1 



c 

0 

1 



c 

0 

1 



10 



globalrow 


1 


1 1 
1 


11 


21 


31 


41 


1 

. 1 . 


51 


61 


71 


81 


1 91 
. 1 


global row 


2 


2 1 
1 


12 


22 


32 


42 


1 

. 1 . 


52 


62 


72 


82 


1 92 
. 1 






-p -p— 1 

3 1 

1 

» • 


13 

4 


23 

» i 


33 

» 4 


43 

1 


1 

-I- 

• 


53 


63 


73 

» 4 


83 


1 93 

.| 


globalrow 


1 

10 1 

1 


1 ^-1 

1 10 1 

1 


4 



20 1 


1 < 

1 

1 30 


» i 

1 1 

1 40 1 


» 

1 

1 50 


• 

•I- 

1 


1 

60 1 


1 1 

1 70 ! 


• i 

1 1 

1 80 


» 

1 

1 90 


.| 

1 100 



Continuing with the example given in Fig. 3.1 and 3.2 we could 
transform the two CSELs into any two adjacent MSELs within the 
master grid. In principle we have only two different choices ; 
either two MSELs in one global row (e.g. : MSEL 1 and 11) or in 
one global column (e.g. : MSEL 1 and 2). 

It turns out that the choice MSEL 1 and 11 results in a lower 
bandwidth of the global linear system (in comparison with the 
choice MSEL 1 and 2). 

Each chosen MSEL corresponds to one isoparametric quadratic 8 
nodal point reference element (see ref .[ 1 ], p . 21/39/104 ) . 

We call a global row rsp. column active, if it exists at least one 
chosen MSEL in that row rsp. column. 
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The reference element numbering is assumed to be as follows : 



Fig. 3.5 : reference element numbering 

8 

jl^ ★ 7 

I I 

2 * * 6 

I I 

3 * * 5 

4 

(3.3) Step 3 j Calculate the 8 ( x , y) -coordinates of each CSEL. 

In the case of the region R given in Fig. 3.1 we get the 
following values : 



Table 3.1 ; 



CSEL (corresponding to MSEL 1) ; 



X y 



nodal point 


1 


0.0 


3.0 




2 


0.0 


2.0 




3 


0.0 


1.0 




4 


0.3827 


0.9239 




5 


0.707 


0.707 




6 


1.8535 


1.8535 




7 


3.0 


3.0 




8 


1.5 


3.0 


CSEL (corre; 


sponding to MSEL 


11 ) : 






X 


y 


nodal point 


1 


3.0 


3.0 




2 


1.8535 


1 .3535 




3 


0.707 


0.707 




4 


0.9239 


0.3827 




5 


1 .0 


0.0 




6 


2.0 


0.0 




7 


3.0 


0.0 




8 


3.0 


1.5 



coordinate 



coordinate 



Because we deal with isoparametric elements, the information given 
in table 3.1 is sufficient to construct the (coarse) mesh shown in 
Fig . 3.2. 
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mesh 



(3.4) Step 4 : Decide how to refine the (coarse) 

defined by the CSELs. 

Finally we have to define the refinement of the (coarse) mesh. Vie 
do this by assigning every active global row rsp. column an 
positive integer value called the number of subrows rsp. 
subcolumns . 

Let us resume the example given by Fig. 3.1 and 3.2. If we would 
define ; 

number of subrows of global row 1 ; 2 

number of subcols of global col 1 ; 2 

number of subcols of global col 2:3, 

the mesh generator would produce the (fine) mesh shown in Fig. 
3.6. 

Fig 3.6 : 




(3.5) Summary of the concept 



Enter during input 



Step 1 : Choose CSELs 

Step 2 : Transform CSELs into MSELs 
Step 3 : Compute coordinates of CSELs 
Step 4 ; Define refinement 



identification numbers 
( X, y) -coordinates 
nr of subrows /subcols 



(4) Documentation 



With the use of nested menus and keyboard filtering the program 
requires no learning on the user's part. The menus are self- 
-explanatory and channel the user in the desired direction. 

Input data is entered from the keyboard corresponding to the last 
chapter . 

(4.1) The main menu 

To run the program one has to load the executable file named 
MESHMENU.EXE into RAM memory by typing MESHMENU followed by a 
carriage return (compare with the source ' meshmenu . c ' , see S.l). 
When the progam is executed, the main menu is displayed : 

Menu 4.1 : 



This is the mesh generator main menu ! 

If you want to HIT 

- create an mesh generator input file 1 

- display an mesh generator input file 2 

- display an mesh generator input file graphically 3 

- change an mesh generator input file 4 

- run the mesh generator (an input file must exist ) 5 

- display an mesh generator output file 6 

- display an mesh generator output file graphically 7 

- change an mesh generator output file 8 

- exit to DOS 9 

Hit a key (1-9) : @ 



The cursor will be located at the character position indicated by 
the ' @ ' symbol . 

As one can see the package consists of three major ports : 

choice 

1. input phase of the (coarse) mesh 1-4 

2. (fine) mesh generation 5 

3. (fine) mesh output 6-7 

Choice 3 is not implemented yet. 

Choice 9 returns control to DOS. 
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(4.2) Input phase of the (coarse) mesh 

Choice 1 : 'create an mesh generator input file ' 

The module INPUTCRT (see 8.2) is executed. 

The user is asked for (compare with chapter 3) 

1. a mesh generator input filename 

2. the total number of CSELs 

3. the identification numbers of the MSELs 

4. the number of subrows rsp. subcolumns of the active global rows 
rsp. global columns 

5. the 8 (x,y) coordinates of each CSEL 

Various explanations and tables are displayed during this input 
process in order to inform the user and verify the data already 
inputted. In addition "correct prompts" are displayed to allow 
changes of incorrect input values. 

When the input phase has finished INPUTCRT 

1. creates an mesh generator input file storing all input data. 

2. writes that file with the inputted filename to the default disk 
drive . 

3. verifies step 2 and 3 by displaying messages on the screen. 

4. returns controlto MESHMENU and the main menu is displayed. 



Choice 2 ; 'display an mesh generator input file' 

The module KEYIMAGE (see 8.3) is executed. 

The user is asked for the filename of the (existing) mesh 
generator input file. 

All input data is displayed : 

- chosen MSELs 

- (x,y) coordinates of corresponding CSELs 

- number of subrows rsp. subcolumns of active global rows rsp. 
global columns . 

Finally control is given to MEStlMENU and the main menu is diplayed. 



Choice 3 ; 'display an mesh generator input file graphically' 

The module KEYGRAF (see 8.4) is executed. 

The user is asked for the filename of the (existing) mesh 
generator input file. 

The coarse mesh is displayed graphically . 

Then control is given to MESHMENU and the main menu is displayed. 



Choice 4 : 'change an mesh generator input file' 

The module KEYCHNGE (see 8.5) is executed. 

The user is asked for the filename of the (existing) mesh 
generator input file. 
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Then the following menu is displayed : 
Menu 4.2 : 



If you want to Hit 

- change coordinates 1 

- change the number of subrows 2 

- change the number of subcolumns 3 

- Exit 4 



If one chooses 1,2 or 3 

- the program channels the user through the desired option. 

- returns control to KEYCHANGE and menu 4.2 is displayed. 

If one chooses 4 

the program returns control to MESHMENU and the main menu is 
displayed . 



(4.2) (fine) mesh generation 
Choice 5 : 'run the mesh generator* 

The module MESHGEN (see 3.6) is executed. 

The user is asked to enter the filename of the existing mesh 
generator input file. 

VThile the fine mesh (connectivity and coordinates) is generated, 
informations are displayed to inform the user about the progress 
of the mesh generation. 

3 files with file extensions ' .con* , ' .coo' and '.loc'are created, 
storing connectivity, coordinates and some additional data of the 
fine mesh. 

The files are written to the default disk drive. 

The program verifies all steps by displaying several messages on 
the screen. 

Finally control is returned to MESHMENU and the main menu is 
displayed . 



(4.3) (fine) mesh output 



Choice 6 ; 'display an mesh generator output file' 



The module MESHIMGE 
The user is asked 
generator input ( 1 ) 
Then connectivity 
control is returned 



(see 8.7) is executed. 

for the filename of the corresponding mesh 

file. 

and coordinates are displayed. Thereafter 
to MESHMENU and the main menu is displayed. 
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Choice 7 : 'display an mesh generator output file graphically' 

The module MESHGRAF (see 8.8) is executed. 

The user is asked for the filename of the corresponding mesh 
generator input (1) file. 

Then the fine mesh is displayed graphically. 

Again control is returned to MESHMENU and the main menu is 
displayed . 



(5) Some remarks on the implemented software techniques 

(5.1) dynamic memory allocation 

Memory is allocated dynamically in all modules - except MESHMENU - 
for all variables if the corresponding storage space exceeds 20 
bytes . 

It is established by employing the C system library routines 

- calloc ( ) and 

- free(). (see ref. [4] p. ) 

(5.2) module communication via files 

To take full advantage of small machines with limited RAM memory, 
information is stored on files. The modules communicate 
exclusively via files. 

Basically every module - except MESHMENU - consists of 3 parts : 

- reading an file 

- manipulating that file 

- creating a (new) file . 

In addition to an economical use of RAM memory, this technique 
allows the user to execute the modules independently. For example 
a user may run the input module INPUTCRT (choice 1), exit to DOS 
(choice 9) and switch off the computer. Because all information is 
savedon a file, he may continue his job at a later time without 
loosing any data that is already inputted. 



(5.3) keyboard filtering 

In order to ensure a high degree of user friendliness all keyboard 
input is filtered. 

Integer and floating point filtering is established by employing 
the (Essential S) C library routines getintO and getfloatO (see 
ref. [5], p.4-46). 

In addition the program checks the filenames inputted by the user; 
- if the user wants to run the module INPUTCRT (choice 1) the 
program asks for an filename with file extension '.key'. If the 
user enters an invalid filename, the module displays a "non 
valid" message and asks for another name. 
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- if the user chooses one of the modules 2, 3, 4, 5, 6, or 7 the 
program asks for an (existing) filename with file extension 
'.key'. On entering a filename with a non valid file extension, 
the routine "behaves" as described above. If the user enters a 
filename with valid file extension, the program verifies whether 
that file does exist on the default disk drive. If the file 
does not exist, the module displays a "does not exist" message 
and asks for another filename. 



(5.4) preparation of screen output in RAM memory 

If the program is executed the user will recognize that the main 
menu in MESHMENU (see menu 4.1) or the messages of INPUTCRT are 
displayed almost instantaneously. 

Menus and messages were edited, displayed once on the screen 
and then saved to the disk. Once they were saved, any program has 
access to these "menu" and "message" files. 

The routine wndfrdsk() (see ref . [ 5 ] , p.4- 1 54) restoressuch a file 
into RAM memory. 

For example the main menu (see menu 4.1) is saved in the file 
window. 0 (compare with (8.1) ). After invoking wndfrdsk() in (8.1) 

the file information is restored into the integer array mess0[]. 
The routine rstwindo( ) (see r e f . [ 5 ] , p . 4 - 1 0 5 ) displays the 
information - now saved in RAM memory - on the text mode screen. 
The screen output information is first copied to the PC's graphics 
RAM memory and then displayed. 

(6) File structures 

(6.1) ' .key ' files 

' .key ' files are 

- created in INPUTCRT 

- read in XEYIMAGE, KEYGRAF, KEYCHNGE, MESHGEN 

- changed in KEYCHNGE 

- formatted files 

structure ; 

total # of MSELs (assume k) 

id. # of MSEL 1-k 

- # of "last" active global row (assume 1) 

- - " - col (assume m) 

- # of subrows of global row 1-1 

- # of subcols of global col 1-m 



x-coord . 


of 


nodal 


point 


1-8 


of 


CSEL 


1 


y-coord. 
• • • 


of 


nodal 


point 


1-8 


of 


CSEl 


1 


” • • • 

x-coord. 


of 


nodal 


point 


1-8 


of 


CSEL 


k 


- y-coord. 


of 


nodal 


point 


1-8 


of 


CSEL 


k 
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(6.2) ' . con ’ files 

' . con ' files are 

- created in MESHGEN 

- read in MESHIMGE, MESHGRAF 

- unformatted files 



structure : 

( comp 
f_m 
l_n_p 

n_p_number 
s elem 



computed 

fine mesh 

local nodal point 

nodal point number in the fine mesh 
(sub) element in the fine mesh ) 



Bytes ; 



2 maximum n_p_n (maxnodenr) 

4 total number of s elem in f m (totelems) 



4 + 
4 + 
4 + 
4 + 
4 + 



2 s_elem number in f_m of 1st comp s elem 
4 id number of MSEL this s_elem belongs to 
6 not yet used 

8 not yet used 

10 n p num of 1 n p 1 of 1st comp s elem 



- 4 + 



24 



n p num of 1 n p 8 of 1st comp s elem 



4 


+ 


• • • 


s elem : 


number 


in 


f m of 


last 


comp s elem 


4 




• • • 


id number of 


MSEL 


this s 


elem 


i belongs to 


4 


+ 


# • • 


not yet 


used 










4 


+ 


• • • 


not yet 


used 










4 




# • • 


n p num 


of 1 


_^_P 


1 of 


last 


comp s elem 


4 


• • 

+ 


24* totelems 


n p num 


of 1 


_"_P 


8 of 


last 


comp s elem 
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(6.3) '.coo' files 



' . coo ' files are 

- created in MESHGEN 

- read in MESHIMGE, MESHGRAF 

- unformatted files 

structure : 



Bytes : 



- 


4 


x-coord . 


of 


n_p 


1 


in 


1st 


comp 


s elem 


- 


8 


y-coord . 


of 


n_p 


1 


in 


1 St 


comp 


s elem 




• • • 

60 


x-coord . 


of 


n_p 


8 


in 


1st 


comp 


s elem 


■ 


64 


y-coord . 


of 


n_p 


8 


in 


1st 


comp 


s elem 




• • ♦ 


x-coord . 


of 


n_p 


1 


in 


last 


comp 


s elem 


- 




y-coord . 


of 


n_p 


1 


in 


last 


comp 


s el em 




• • • 


x-coord . 


of 


n_p 


8 


-in 


last 


comp 


s elem 


- 


totelems * 64 


y-coord . 


of 


n_p 


8 


in 


last 


comp 


s elem 
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(6.4) * . loc ' files 



' . loc ' files are 

- created in MESHGEN 

- read in MESHIMGE, MESHGRAF 

- unformatted files 

structure : 

{ loc = location ) 

Bytes ; 



2 


loc 


of 


s elem 


1 


in 


' . con ' / ' • coo ' 


file 


* • « • 

- totelems * 2 


loc 


of 


s elem 


totelems 


in 


' . con ' / ' • coo ' 


file 
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(8) Source Listings 

(8.1) meshmenu.c 

finclude <stdio.h> 
#define TRUE 1 



main ( ) 

{ 

int mess0[ 1541 ] , j ; 

wndfrdsk ( "window . 0 ", mess0 ) ; /* meshmenu main display */ 

do { 

clrscrn ( ) ; 
clearkbd ( ) ; 

rstwindo (l,5,23,71,mess0) ; 

curlocat (22,24) ; 

get int (l,l,0,l,2,&j,l,l,9); 

curlocat ( 25 , 0 ) ; /* cursor vanishes */ 

clrscrn ( ) ; 



switch(j) { 



case 1 

case 2 

case 3 

case 4 

case 5 

case 6 

case 7 

case 8 

case 9 
default 

} /* end 



: inputcrt ( ) ; 
break; 

: keyimage ( ) ; 
break ; 

; keygraf ( ) ; 
break ; 

: keychnge(); 
break ; 

: meshgen ( ) ; 
break ; 

; meshimge ( ) ; 
break ; 

; neshgraf ( ) ; 
break; 

; printf ( " ' meshchange . c ' does not exist yet 
break ; 

; return; 

; printf ("You pressed a non valid key 1“); 
break ; 

switch */ 



1 ") 



if ( ( j 1=2 )&& ( j 1=3 )&& ( j 1=4 )&& ( j 1=6 )&&( j ! =7 ) ) { 

curlocat (24,0) ; 

Printf ("Hit any key to continue ..."); 
getkey (S( j ) ; 



} while (TRUE) ; 



} 
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(8.2) inputcrt.c 



tinclude <stdio.h> 
#include <malloc.h> 
static char s[15]; 
static char *name=s; 



inputcrt ( ) 

{ 

FILE * input; 

unsigned short *keygrid, *idnumber, *nr_rows , *nr_cols, lastcol=0, 
lastrow=0 , i=0 , k, row, col , nr, rowml , colml , nr_key elements ; 
float *x CO key , *y_co_key ; 

int *messl , *mess2 , *mess3 , mess4[ 2] , messblkl [l],messblk2[2]; 



messl = (int * ) calloc ( 133 , sizeof ( int ) ) ; 
mess2 = (int *) calloc ( 396 , sizeof ( int )) ; 
mess3 = (int *) calloc ( 1584 , sizeof ( int )) ; 

keygrid = (unsigned short *) calloc ( 101 , sizeof ( unsigned 
idnumber = (unsigned short *) calloc ( 101 , sizeof ( unsigned 
nr_rows = (unsigned short *) calloc (11 , sizeof ( unsigned 
nr cols = (unsigned short *)calloc(ll , sizeof ( unsigned 



short ) ) ; 
short ) ) ; 
short ) ) ; 
short ) ) ; 



if 



( raessl==NULL 1 1 
idnumber==NULL I 
printf ("*calloc() 
return; 

1 



mess2==NULL 1 1 mess3==NULL 1 I keygrid==NULL 1 1 
1 nr_rows==NULL I 1 nr_cols==NULL ) { 
fails, not enough storage space I"); 



wndf rdsk ( "wndblk . 1 , messblkl ) ; 
wndf rdsk ( "wndblk . 2 " , messblk2 ) ; 
wndf rdsk ( "window. 1 " , messl ) ; 
wndf rdsk ( "window. 2 " , mess2 ) ; 
wndf rdsk ( "window. 3 " , mess3 ) ; 
wndf rdsk ( "window. 4 " , mess4 ) ; 



/* 


* 


messblkC 1 ] 


*/ 


/* 


** 


messblkC 2] 


*/ 


/* 


node connec . 


messl [ 133 ] 


*/ 


/* 


small keygrid 


mess2[396] 


*/ 


/* 


large keygrid 


mess3 C 1 584 ] 


*/ 


/* 


2 blanks 


mess4[ 2 ] 


*/ 



clrscrn ( ) ; 

printf ( "During this input procedure the program creates an\n"); 
printf (" input file l\n"); 

printf ( "\nEnter the name of the inputfile to be created\n"); 
printf ("( Use the file extension '.key' ) : "); 

label; scanf ("%s",s); 



wh ile(s[i]l='\0'){ 

if (s[i] 1= '.' ) { i=i+l ; continue ; } 

break; 

] I 

if((s[i+l] i= 'k' ) I I (s[i+2]! = 'e' ) 1 I (s[i+3]l = 'y' ) ) { 

printf ("The inputfile has a non valid file extension '.key'\n")| 
printf ("Try again : "); ; 

goto label; 



input = fopen ( name , "w" ) ; 
clrscrn( ) ; 

rstwindo(0, 0,23, 6 5, mess 3 ) ; 
curlocat (24,0) ; 

printf("Hit any key to continue 
getkey ( &i ) ; 

do { 

clrscrn ( ) ; 

printf ( "\nEnter the number of keyelements : "); 

get int (3, 1,-1, 1,2, &nr_keyelements ,0,1,100); 

curlocat (24,0); 

printf ( "correct ? (Y,N) "); 

i = ecoyesno ( 24 , 16 , 1 ) ; 

} while ( i==0) ; 

fprintf ( input , " %u " , nr_keyelements ) ; 

x_co_key = (float *) calloc ( 8*nr_keyelements , si zeof ( float )) ; 
y_co_key = (float *) cal loc ( 8*nr_keyelements , sizeof ( float )) ; 

if (x_co_key==NULL I I y_co_key==NULL ) { 

printf ( "*calloc ( ) fails, not enough storage space 1"); 
return ; 

} 

clrscrn ( ) ; 

rstwindo (1,47, 12,79,mess2) ; 
for ( nr=l ; nr <=nr_keyelements ; ++nr ) { 

do { 

curlocat (15,0); 
clreos ( ) ; 
curlocat (15,0); 

print f ( "Enter id. number of (first or) next keyelement 

ge tint (2, 1,-1, 1,2, &idnumber[ nr ] ,0,1,99) ; 

curlocat (24,0); 

printf ( "correct ? (Y,N) "); 

i = ecoyesno ( 24, 16 , 1 ) ; 

} while ( i == 0 ) ; 

fprintf ( input , " %u " , idnumber[nr ] ) ; 

keygridC idnumber [ nr ] ] = 1; 

rstwindo ( rowml + 1 , 45 , rowinl + 1 , 40 , mess4 ) ; 

rstwindo ( 0 , 46 + 3 *colml , 0 , 47 + 3*colml , mess4 ) ; 

col = idnumber[nr]/l0 . 5 + 1; 

row = idnumberCnr] - (col-l)*10; 

paint ( row+1 , 46+3*col , row+1 , 48+3*col ,23,0); 

rstwindo ( row+1 , 45 , row+1 , 45 , messblkl ) ; 

rstwindo (0, 46+3* col, 0,47+3* col, messblk2); 

colml =- col ; 

rowml = row; 

if ( lastrow < row ) lastrow = row; 
if ( lastcol < col ) lastcol = col; 

) /* end nr */ 
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curlocat ( 25 , 0 ) ; /* cursor vanishes */ 

rstwindo{ rowml+1 , 45 , rowml+1 , 46 ,mess4 ) ; 
rstwindo ( 0 , 46+3*colral , 0 , 47+3*colml , mess 4 ) ; 

savwindo{ 1 , 47 , 12 , 79 ,raess2 ) 7 /* saving the blinking grid */ 
fprintf ( input , " %u ",lastrow); 
fprintf ( input , " %u ", lasted); 

for ( row=l 7 row<=lastrow; ++row) { 
curlocat ( row+1 , 45 ) ; 
putchar ( 16 ) ; 

paint ( row+1 , 45 , row+1 , 45 , 23 , 0 ) 7 /* makes it blink */ 

do { 

curlocat ( 12+row, 0 ) ; 
clreos ( ) 7 

curlocat ( 12+row, 0 ) ; 

print f ( "Enter number of subrows in globalrow %u : ",row); 

clearkbd { ) ; 

getint ( 2 , 1 , -1 , 1 , 2 , &nr_rows[ row] ,0,1,99); 
curlocat (24,0); 
printf( "correct ? (Y,N) "); 

i = ecoyesno { 24 , 16 , 1 ) 7 
} while (i==0); 

fprintf ( input , " %u " ,nr_rows[row] ) ; 
rstwindo( row+1 , 45 , row+1 , 46 , mess4 ) ; 

) 

for { col=l 7 col <=lastcol 7 ++col ) { 

curlocat ( 0 , 46+3*col ) ; 
putchar ( ' \31 ' ) ; 

paint (0,46+ 3* col, 0,46+3* CO 1,23,0) ; 
do { 

curlocat ( 1 2+col , 0 ) ; 
clreos ( ) 7 

curlocat ( 1 2+col , 0 ) ; 

printf ( "Enter number of subcolumns in globalcolumn %u : ",col) 

getint (2, 1,-1, 1,2, &nr_cols[col] ,0,1,99); 
curlocat { 24 , G ) ; 
printf ( "correct ? (Y,N) "); 

i = ecoyesno ( 24 , 16 , 1 ) 7 
} while ( i==0 ) 7 

fpr intf ( input , " %u " , nr_cols[ col ] ) ; 
rstwindo ( 0 , 4 6+ 3* col , 0 , 47+3* col , mess4 ) ; 

} 

clrscrn ( ) ; 

printf ("The system will now ask for coordinates l\n\n\n"); 
printf ("The node connec . within the (reference) keyelements\n" ) ; 
printf ("is as follows :"); 
rstwindo ( 6 , 13,12,31,messl) ; 
curlocat ( 24 , 0 ) ; 

printf ("Hit any key to continue ..."); 
syspause ( 0 , 0 , 1 , 0 ) ; 
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getkey ( &i ) ; 
clrscrn( ) ; 

rstwindo(0, 60,6, 78,messl ) ; 

print f ( "Enter the node coordinate at actual cursorposi tion \n"); 
printfC'and hit 'ENTER'. " ) ; 
curlocat (7,0); 

for ( nr=l ; nr <=nr_keyelements ; ++nr ) { 

col = idnumber[nr]/l0 . 5 + 1; 
row = idnumberCnr] - (col-1 )*10; 



printf ( " identification 


number 


; %3u 


\n " , idnumberCnr] ) ; 


printf (" in 


global 


column 


: %3u 


\n" ,col ) ; 


printf ( "and 


global 


row 


: %3u 


\n\n" , row) ; 


printf ( " 


x-coordinate 


y- 


■coordinate \n"); 



for ( k=l ; k<=8 ; ++k ) { 

do { 

curlocat ( 11+k, 0 ) ; 
clreos ( ) ; 

rstwindo (8,47,19,79, mess 2 ) ; 

rstwindo ( row+8 , 45 , row+8 , 45 , raessblkl ) ; 

rstwindo (7,46+3*col,7,47+3*col,messblk2); 

curlocat ( 11+k, 0 ) ; 

printf("node %lu ; ",k); 

getf loat(6,6, 1, -1, 1 , 2 , &x_co_key[ ( nr-1 ) *8+k] , 0 , -1008 . , 1000 . ) 
curlocat ( 1 1+k , 27 ) ; 

getfloat(6,6, 1,-1, 1,2, &y_co_key[ ( nr-1 ) *8+k] , 0 , - 1O0O . , lOOO . ) 

curlocat (24,0); 

printf ( "correct ? (Y,N) "); 

i = ecoyesno( 24, 16 , 1 ) ; 

} while ( i==0 ) ; 

} /* end k */ 

for (k=l ;k<=8 ; ++k) fprintf ( input , " %f " , x_co_key[ ( nr-1 ) *S+k] ) ; 
for (k=l ; k<=8 ; ++k ) fpr intf ( input , " %f " , y_co_key[ ( nr-1 ) *S+k] ) ; 
curlocat (7,0); 
clreos ( ) ; 

) /* end nr */ 
fclose ( input ) ; 
f ree (messl ) ; 
free (mess 2 ) ; 
free ( mess3 ) ; 
f ree ( keygrid ) ; 
free ( idnumber ) ; 
free ( nr_rows ) ; 
f ree ( nr_cols ) ; 
free ( x_co_key ) ; 
f ree ( y_co_key ) ; 
clrscrn( ) ; 

printf ("The inputfile with name ; %s \n",s); 

printf ("is created and saved on the actual discdr ive . \n\n\n " ) ; 

1 /* end inputcrt */ 
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(8.3) keyimage.c 



#include <stdio.h> 

#include <malloc.h> 
static char s[15]; 
static char *name = s; 

keyimage ( ) 

{ 

FILE *input; 

unsigned short *keygrid, *idnumber, lasted , lastrow, j ,k, row, col , 

*nr_rows , *nr_cols , nr, nr_keyelements ; 
float *x_co_key , *y_co_key ; 
int *messl , *mess2 ; 

messl = (int * ) calloc ( 133 , sizeof ( int ) ) ; 
mess2 = (int * )calloc( 396, sizeof ( int )) ; 

keygrid = (unsigned short *) calloc ( 101 , sizeof ( unsigned short)); 
idnumber = (unsigned short *) calloc ( 101 , sizeof ( unsigned short)); 
nr_rows =( unsigned short *) calloc ( 11 , sizeof ( unsigned short)); 
nr cols = (unsigned short *) calloc ( 11 , sizeof ( unsigned short)); 



clrscrn ( ) ; 

if ( messl==NULL I I mess2==NULL I I keygrid==NULL I I 

idnumber==NULL I I nr rows==NULL I I nr cols==NULL ) 

{ 

printf(" *calloc() fails, not enough storage space I"); 
return; 

} 

wndfrdsk ( "window. 1 ", messl ) ; /* node connectivity */ 

wndfrdsk( "window. 2" ,mess2 ) ; /* small keygrid */ 

printf("This subroutine displays the data of an existing \n 
printf ( "meshgenerator inputf ile . \n\n" ) ; 
print f ( "Enter name of inputf ile : "); 

scanf ( "%s " , s ) ; 

v/hile ( f ilexist ( name ) == 0 ) { 

printf (" \nThe file %s does not exist I Try again ; ",s); 
scanf ( "%s " , s ) ; 

} 

input = f open ( name ," r ") ; 
f scanf ( input, "%u" , &nr_keyelements ) ; 

x_co_key = (float *) calloc ( 8*nr_keyelements , sizeof ( float )) ; 
y_co_key = (float *) calloc ( 8*nr_keyelements , sizeof ( float )) ; 

if ( x_co key==NULL 1 I y_co_key==NULL ) { 

printfT"*calloc ( ) fails for x_co_key/y_co_key 1"); 
return ; 

} 
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for ( nr=l ; nr <=nr_keyelements ; ++nr ) { 

fscanf ( input, "%u" , & idnumber [ nr ] ) ; 
keygridC idnumber[nr ] ] = 1; 

} 

fscanf { input, "%u" , Sclastrow) ; 
fscanf { input, "%u" , &lastcol ) ; 

for ( row=l ; row<=lastrow; ++row) fscanf ( input , " %u " , & nr rows [ row] ) 
for ( col=l ;coK=lastcol; ++col ) f scanf ( input , "%u " , &nr_col s[ col ] ) 
for ( nr=l ; nr <=nr_keyelements ; ++nr ) { 

for (k=l ; k <=8 ; ++k ) f scanf ( input ,"% f" , &x_co_key[ ( nr-1 ) *8+k ] ) 
for ( k=l ; k <=8 ; ++k ) fscanf ( input , "%f" , &y_co_key[ ( nr-1 ) *8+k] ) 
] /* end nr */ 
fclose ( input ) ; 



/* image of input */ 



clrscrn ( ) ; 

printf("Here the display of 
printf ( "meshgenerator input 
printf ( " %s " , s ) ; 

curlocat (9,4) ; 
printf ("( active keyelements 
curlocat (10,4); 
printf("(are blinking. 



the\n" ) ; 
file : \n\n\n" 

) ■• ) ; 

) " ) ; 



rstwindo( 1,43,12,75, mess 2 ) ; 
for ( nr=l ; nr <=nr_keyelements ; ++nr ) { 

col = idnumber [ nr ] /10 . 5 + 1; 
row = idnumber [nr] - (col-1 )* 10; 
paint ( row+1 , 42+3*col , row+1 , 44+3*col ,23,0) ; 
] 



curlocat (0,64) ; 

printf ("nr of subrows "); 

putchar ( ' \31 ' ) ; 

for ( row=l ; row< = lastrow; ++rov/) { 

curlocat ( row+1 ,77); 
print f (" %2u ", nr rows[row]); 

} 

curlocat(14,46) ; 

putchar ( ' \30 ' ) ; 

printf(" nr of subcoluinns"); 

curlocat (13,44) ; 

for ( col = l ; coK = lastcol ; ++col ) printf ( "°3u" , nr_cols[col ]) ; 

for ( nr=l ; nr <=nr_keyelements ; ++nr ) { 

col = idnumber [ nr ]/ 10 . 5 +1; 
row = idnumber[nr] - (col-1) *10; 
curlocat(15,0) ; 

rstwindo (16,60,22,78, mess 1 ) ; 

printf (" idnumber of keyelement : ° 2u\n" , idnumber[nr] ) ; 
for (k=l;k<=8;++k) printf("x%lu : %12f y%lu : ^12f\n", 

k,x CO keyC ( nr-1 ) *8+k] , k , y_co_keyC ( nr- 1 ) *0+k] ) ; 
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curlocat (24,50); 

printf{"Hit any key to continue 
getkey (& j ) ; 
curlocat (15,0); 
clreos ( ) ; 

} /* end nr */ 

free ( messl ) ; 
free (mess2 ) ; 
free(keygrid) ; 
free ( idnumber ) ; 
f ree ( nr_rows ) ; 
free ( nr_cols ) ; 
free ( x_co_key ) ; 
free ( y_co__key ) ; 

} /* end keyimage */ 
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(8.4) keygraf.c 



#include <stdio.h> 

#include <malloc.h> 
char s[15]; 
char *name=s; 

keygraf ( ) 

{ 

FILE *input; 

unsigned short nr_keyelements , nr , j , k, factor_x , lastrow, lasted , 

nr_rows , nr_cols , row, col, idnumber, factor y, factor; 
int *x_int, *y_int ; 

float *x, *y, *x_co_key, *y_co_key, *x key , *y key; 

float min_x= 10000 . , min_y= 100000 . , max_x=- 10000 . , max_y=- 10000 . ; 

x_int = (int * ) calloc ( 25 , sizeof ( int ) ) ; 
y_int = (int *) calloc ( 25 , sizeof ( int )) ; 

X = (float *) calloc ( 25 , sizeof ( float )) ; 
y = (float *) calloc ( 25 , sizeof ( float )) ; 
x_key = (float *) calloc ( 9 , sizeof ( float )) ; 
y_key = (float *) calloc ( 9 , sizeof ( float )) ; 

if (x_int==NULL I I y_int==NULL I I x==NULL I | y==NULL I I 
X key==NULL I I y key==NULL ) 

{ “ 

print f (" *calloc ( ) fails 1"); 
return; 

) 



clrscrn ( ) ; 

print f ( "Enter the input filename to be displayed graphically : "); 

scanf ( "%s" , s ) ; 

while ( f ilex ist ( name ) == 0 ) ( 

printf ( "\nThe file %s does not exist 1 Try again : ",s); 

scanf ( "%s" , s ) ; 

} 



/* read inputdataf ile */ 

input = fopen( name , " r " ) ; 

fscanf ( input, "%u" , Stnr_keyelements ) ; 

for ( nr=l ; nr <=nr_keyelements ; ++nr ) f scanf ( input %u " , & idnumber ) 
fscanf ( input, "%u" , Sdastrow) ; 
fscanf ( input, "%u" , Sclastcol ) ; 

for ( row=l ; row <=las trow; ++row) f scanf ( input , " ° u " , &nr_rows ) ; 
for ( col=l ; col <=lastcol ; ++col ) f scanf ( input , "%u " , &nr_col s ) ; 

j=8*nr_keyelements +1; 

x_co_key =( float *) calloc ( j , sizeof ( float )) ; 
y CO key =(float *) calloc ( j , sizeof ( float )) ; 
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if ( X CO key==NULL I 1 y co key==NULL ) 

printf ( "*calloc ( ) fails I"); 
return ; 

} 



for ( nr=l ; nr <=nr_keyelements ; ++nr ) { 

for (k=l ; k<=8 ; ++k ) f scanf ( input , " %f " , Scx_co_key[ ( nr-1 ) *8+k ] ) ; 
for (k=l ; k<=8 ; ++k ) f scanf ( input ,"% f" , &y co key[ ( nr-1 ) *8+k] ) ; 
} " " 
fclose ( input ) ; 

for (k=l ; k<=8*nr_keyelements ; ++k ) { 

if (min_x > = x_co_key[k] ) min_x = x_co_key[k]; 

if (min_y >= y_co_key[k] ) min_y = y_co_key[k]; 

if (max_x < = x_co_key[k] ) max_x = x_co_key[k]; 

if (max y <= y CO key[k] ) max y = y co key[k]; 

} - - - _ _ _ 
factor_x = 303 . 0/ (max_x-min_x ) ; 
factor_y = 183 .0/ (max_y-min_y ) ; 
if ( factor_x >= factor_y ) factor = factor_y; 
else factor = factor_x; 

grafinit(2,0,0) ; 

for ( nr=l ; nr <=nr_keyelements ; ++nr ) ( 

for ( j=l ; j <=8; ++ j ) { 

x_key[j] = x_co_key[ ( nr-1 ) *8+ j ] ; 
y keyCj] = y co key[ ( nr-1 ) *8+ j ] ; 

}_ _ _ 
for ( j=l ; j<=8;++j ) { 

x[l+(j-l)*3] = x_keyCj]; 
y[l+(j-l)*3] = y key[j]; 

} 

for ( j=l; j<=2;++j) { 

quad_2d(-1.0 ,1.0 - j /3 . 0 , x_key , y_key , &x [ j+ l],&yCj+ 1]); 

quad_2d(-1.0 , - ( j/3 . 0 ) , x_key , y_key , &x[ j+ 4],&y[j+ 4]); 

quad_2d ( -1 . 0+ j/ 3 . 0 , -1 . 0 , x_key , y_key , &xC j+ 7],&y[j+ 7]); 

quad_2d( j/3.0,-1.0 , x_key , y_key , &x[ j+10 ] , &y[ j+13 ] ) ; 

quad_2d( 1.0 , -1 . 0+ j/ 3 . 0 , x_key , y_key , &x[ j + 13 ] , £<y[ j + 1 3 ] ) ; 

quad_2d( 1.0 , j / 3 . 0 , x_key , y_key , &x[ j+16 ] , &y [ j+16 ] ) ; 

quad_2d( 1.0-j/3.0, 1.0 , x_key,y_key, &x[ j+19] , &y[ j+19] ) ; 

quad 2d( -(j/3.0), 1.0 ,x key,y key , &x[ j+22] , &y[ j+22 ] ) ; 

}_ _ _ 
for ( j=0; j <=23;++ j ) { 

x_int[j] =3 + factor * ( x[j+l] - min_x ); 

y intCj] = 191 - factor * ( y[j+l] - min y ); 

}“ 

grpoly( x_int, y_int , 24,1); 

} /* end nr */ 

grxlab("Hit any key to continue 7 , 55 , 0 ) ; 

getkeyUj); 

graf init (0,0,0); 



I 
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free { x_int ) ; 
free { y int ) ; 
free ( xT; 
free(y ) ; 
free { x_key ) ; 
free(y_key) ; 
free { x_co_key ) ; 
free ( y_co_key ) ; 

} 
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(8.5) keychnge.c 



#include <stdio.h> 
finclude <malloc.h> 
static char s[15]; 
static char *name = s; 

keychnge ( ) 

{ 

FILE * input; 

unsigned short *keygrid ,* idnuraber , lasted , lastrow, j , k, row, col , 

*nr_rows, *nr_cols, idnr , nr, nr_keyelements ; 
float *x_co_key , *y_co_key ; 
int i , *messl , *mess2 ; 

messl = (int * ) calloc ( 133 , sizeof ( int ) ) ; 
mess2 = (int *) calloc ( 396 , sizeof ( int )) ; 

keygrid = (unsigned short *) calloc ( 101 , sizeof ( unsigned short)); 
idnumber = (unsigned short *) calloc ( 101 , sizeof ( unsigned short)); 
nr_rows =(unsigned short *) cal loc ( 11 , sizeof ( unsigned short)); 
nr cols = (unsigned short *) calloc ( 11 , sizeof ( unsigned short)); 



clrscrn ( ) ; 

if ( raessl==MULL 1 I mess2==NULL I | keygrid==NULL I I 

idnumber==NULL I I nr rows==NULL I I nr cols==NULL ) 

{ 

printf(" *calloc() fails, not enough storage space 1"); 
return ; 

} 

wndfrdsk ( "window. 1 ", messl ) ; /* node connectivity */ 

wndfrdsk ( "window. 2 ", mess2 ) ; /* small keygrid */ 

printf("This subroutine allows you to change an existing\n") 
printf ( "meshgenerator inputf ile . \n\n" ) ; 
printf ( "Enter name of inputfile : "); 

scanf ( "%s " , s ) ; 

while ( f ilexist ( name ) == 0 ) { 

printf ( "\nThe file %s does not exist 1 Try again : ",s) 

scanf ( "%s " , s ) ; 

} 

input = f open ( name, "r" ) ; 
f scanf ( input, "%u" , &nr_keyelements ) ; 

x_co_key = (float *) calloc ( 8*nr_keyelements+l , sizeof ( float )) ; 
y_co_key = (float *) calloc ( 8*nr_keyelements+l , sizeof ( float )) ; 

if ( x_co_key==NULL I I y_co_key==NULL ) { 

print f (" *cal loc ( ) fails for x_co_key/y_co_key 1"); 
return; 

} 
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for ( nr=l ; nr <=nr_keyelements ; ++nr ) { 

fscanf ( input, "%u“ , &idnumber[nr] ) ; 
keygrid[idnumber[nr ] ] = 1; 

} 

f scanf ( input , "%u " , & las trow) ; 
fscanf ( input, "%u" , &1 as tool ) ; 

for ( row=l ; row< = lastrow; ++row) fscanf ( input , " %u ” , &nr_rows [ rov/] ) 
for ( col=l ;coK=lastcol; ++col ) f scanf { input , " %u" , &nr_cols[ col ] ) 
for ( nr=l ; nr <=nr_keyelements ; ++nr ) { ~ 

for (k=l ; k<=8 ; ++k ) fscanf ( input , "%f" , &x_co_key[ { nr-1 ) *8+k] ) 
for (k=l ; k<=8 ; ++k ) f scanf { input ," %f" , &y_co_key[ ( nr-1 ) *8+k ] ) 
} /* end nr */ 
f close ( input ) ; 



do { 

clrscrn ( ) ; 

printf(" If you want to 
printf(" change coordinates 
printf(" change the number of subrows 
printf(" change the number of subcolumns 
printf(" Exit 

printf ( " \nHit a key (1-4) : "); 

curlocat(7, 18) ; 
getint( 1 ,l,0,l,2,&j,l,l,8) ; 
clrscrn ( ) ; 



HIT \n\n" ) ; 

1 \n"); 

2 \n"); 

3 \n"); 

4 \n“); 



switch( j ) { 

case 1 : clrscrn(); 

rstwindo (1,47, 12, 79, me ss2) ; 

for ( nr=l ; nr <=nr_keyelements ; ++nr ) 

{ /* keygrid blinks */ 

col = idnumber[nr]/l0 . 5 + 1; 

row = idnumberCnr] - (col-1 )*10; 

paint ( row+1 , 46 + 3 * col , row+1 , 48+3 *col , 23 , 8 ) ; 

) 

do { 

curlocat(14,0) ; 
clreos ( ) ; 
curlocat(14,0) ; 

printfC'Enter id# of element to be changed ;"); 
curlocat(14,41) ; 

getint (2, 1,-1, 1,2, 5<idnr , 0,1,99); 
curlocat (24,0) ; 
printf ( "correct ? (Y,M) "); 

i = ecoyesno ( 24 , 16 , 1 ) ; 

} while (i == 0); 



for (nr=l;nr<=nr keyelements ; ++nr ) 
if ( idnumberrnr ]==idnr ) break; 
clrscrn( ' ; 

printf(" old new "); 

printf (" correct \n\n"); 

for (k=l ; k <=8 ; ++k ) 
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printf("x%lu : %10f y%lu : %10f\n" , k, x_co_key[ ( nr-1 ) *8+k] 

,k,y_co_key[ (nr-1 )*8+k,k] ) ; 

for (k=l ;k<=8;++k) { 
do { 

curlocat ( k+1 , 69 ) ; 

printf ( "%1 u(Y/nW : ",k); 

i=ecoyesno(k+l ,79,1); 

if (i==0) { 

cur locat { k+1 , 38 ) ; 

clreol ( ) ; 

curlocat (22,0); 

printf ( "Enter new coord, at cursorposition\nand hit 'ENTER'"); 
curlocat (k+1 , 34 ) ; 
printf("x%lu : ",k); 
curlocat ( k+1 , 39 ) ; 

getfloat(6,6, 1,-1, 1, 2,&x_co_key[ (nr-1 )*8+k] ,0,-1000. , 1000 . ) ; 
curlocat (k+1 , 54 ) ; 
printf ("y%lu : ",k); 
curlocat (k+1 , 59 ) ; 

getfloat(6,6, 1,-1, 1, 2,&y_co_key[ ( nr-1 ) *8+k] , 0 , -1000 . , 1000 . ) ; 
curlocat (22,0); 
clreos ( ) ; 

} /* end if */ 

} while(i==0); 

] /* end k */ 
break; 



case 2 ; 
clrscrn ( ) ; 

printf (" old new "); 

printf (" correct\n\n" ) ; 

for ( row=l ; row<=lastrow; ++row) 

printf ("nr of subrows in row %lu ; %5u\n" , row, nr_rows[ row] ) ; 
for ( row=l ; row<=lastrow; ++row) { 
do { 

curlocat ( row+1 , 54 ) ; 
printf (" (Y/ n)? ; "); 

i=ecoyesno ( row+1 ,63,1); 
if (i==0) { 
curlocat ( row+1 , 36 ) ; 
clreol ( ) ; 
curlocat (22,0); 

printf ( "Enter new 'nr of subrows' at actual cursorposition" ) ; 
printf ( "\nand hit 'ENTER'"); 
curlocat ( row+1 , 36 ) ; 

getint (2, 1,-1, 1,2, &nr_rows[ row] ,0,1,99); 
curlocat ( 22 , 0 ) ; 
clreos ( ) ; 

} /* end if */, 

} wliile ( i==0 ) ; 

} /* end row */ 
break ; 
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1 



old 



new 



case 3 : 
clrscrn( ) ; 
printf ( ’’ 

printf(" correct\n\n" ) ; 

for ( col=l ; col <=lastcol ; ++C0I ) 
printf (“nr of subcols in col %lu : 
for (col=l;coK=lastcol;++col) { 

do { 

curlocat ( col+1 , 54 



%5u\n" , col , nr 






colsCcol ] ) 



printf C’(Y/N)? : "); 

i=ecoyesno ( col+1 ,63,1); 
if (i==0) { 

curlocat ( col + 1 , 36 ) ; 
clreol ( ) ; 
curlocat (22,0) ; 

printf ("Enter new 'nr of subcols’ at actual cursorposition" ) ; 
printf (’’ \nand hit 'ENTER'"); 

curlocat(col+l, 36) ; 

getint(2, 1,-1, 1,2, Scnr_cols[col] ,0,1,99); 
curlocat ( 22 , 0 ) ; 
clreos ( ) ; 



} 

} while ( i==0 ) ; 

) 



break ; 

case 4 : return; 
default ; break; 

} /* end switch */ 

input = f open ( name , "w" ) ; 
fprintf (input, " %u " , nr_keyelements ) ; 
for ( nr=l ; nr <=nr_keyelements ; ++nr ) 

fpr intf ( input , " %u " , idnumberC nr ] ) ; 
fprintf ( input , " %u ",lastrow); 
fprintf Unput , " %u ",lastcol); 
for ( row=l ; row<=lastrow; ++row) 

fprintf ( input, " %u " , nr_rows[ row] ) ; 
for ( col=l ; col <=lastcol ; ++C0I ) 

fprintf ( input, " %u " , nr_cols[col] ) ; 
for ( nr=l ; nr<=nr_keyelements ; ++nr ) { 

for (k=l ;k<=8;++k) 

fprintf ( input , " ?,f " , x_co_key[ ( nr-1 ) *8+k] ) ; 

for (k=l ;k<=8;++k) 

fprintf ( input , " %f " , y_co_key[ ( nr- 1 ) *8+k] ) ; 

} 

fclose ( input ) ; 

1 while (1); /* global do */ 

f ree(messl ) ; 

free (mess2 ) ; 

free ( keygrid) ; 

free ( idnumber ) ; 

free ( nr_rows ) ; 

free ( nr_cols ) ; 

free ( x_co_key ) ; 

free(y_co_key) ; 

) /* end keychnge */ 
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(8.6) meshgen.c 

#include <stdio.h> 

#include <raalloc.h> 

static char s_loc[ 15 ] , s_conn[ 15 ] , s_coord[ 15 ] ; 

static char *nameloc=s_loc , *nameconn = s_conn , *namecoord=s_coord ; 
static float *coord; 

static unsigned short * loc , *connec ; 

meshgen ( ) 

{ 

FILE *locfile, *connfile, *coordfile; 

unsigned short las trow, las t col , nr_of_sube lament s[ 11], nr_rows[ 11], 

nr_cols[ 11 ] , *keygrid , totelems , maxnodenrs=0 ; 
unsigned short glcol_f irstnode=l , glcol_subelement=l , glcol , glrow, 

key_f irstnode , key_subelement , f irstnode , subelement , 
subcol_f irstnode, subcol_subelement , subcol , subrow, 
elemjump, key_col jumpl , key_col jump2 , col jumpl , 
col jump2 , memoryl , memory2 , helpl , help2 , 
noconnect , noconnl , noconn2 , row, col , j , k, 
nr_of_keyelement=0 , gl_nodenr_of_loc_node[9] , 
count 1=1 , count2=0 , dummy=0 , idnumber ; 
float *x_co_key,*y co_key, x[9] , y[9] ,ksi[4] , eta[4] , 

x_co_keyelemr9 ] , y_co_keyelem[ 9 ] ; 

(unsigned short *) calloc ( 1201 , sizeof ( unsigned short)); 

( float * ) calloc ( 1601 , sizeof (float ) ) ; 

(unsigned short *) calloc ( 101 , sizeof ( unsigned short)); 

= (float *) calloc ( 400 , sizeof ( float )) ; 

= (float *) calloc ( 400 , sizeof ( float )) ; 



meshinputread ( s_loc , s_conn, s_coord , x_co_key , y_co_key , keygr id , 

nr_rows , nr_cols , Sctotelems , nr_of_subelements , &lastrow,,Sclastcol ) ; 

loc =(unsigned short *) cal loc ( totelems+1 , sizeof ( unsigned short)); 

if ( connec==NULL 1 1 coord==NULL I 1 loc==NULL 1 I 

keygrid==NULL II x co key==NULL I I y co key==MULL ) 

{ ” “ " ■ 
printf ( "*calloc ( ) fails in meshgen 1"); 
return ; 

} 



connfile = fopen ( naraeconn , "w+b" ) ; 

coordfile = f open ( namecoord, "w+b" ) ; | 

fwrite (( unsigned short *) &maxnodenrs , sizeof ( unsigned short ), 1 , conn f;l 
fwrite (( unsigned short *)&totelems , sizeof ( unsigned short ) , 1 , connf I 

i 

/* Begin of actual computation */ 

for ( glcol=l ; glcol <=lastcol ; ++glcol ) { 

elemjump = 0; 

coljumpl = 0; I 






connec = 
coord = 
keygrid = 
x_co_key 
y CO key 
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coljump2 = 0; 
noconnect = 1 ; 

for ( row=l ; row<=lastrow; ++row) 

if (keygridC (glcol-1 ) *10+row]==l ) { 

elemjump = elemjump + nr_rows[ row] ; 

coljumpl = coljumpl + 2*nr_rows[ row] + noconnect; 

coljump2 = coljump2 + nr_rows[row] + noconnect; 

noconnect = 0; 

) 

else noconnect = 1; 
key_f irstnode = glcol_f irstnode; 
key subelement = glcol subelement; 
key_coljumpl = coljumpl; 
key_coljump2 = coljump2; 

for ( glrow=l ; glrow<=lastrow; ++glrow) { 
idnumber = (glcol-1 )*10 + glrow; 
if ( keygridC idnumber ] == 1 ) { 

nr_of_k eye 1 erne nt = nr_of_keyelement +1; 

subcol_f irstnode = key_f irstnode ; 
subcol subelement = key_subelement ; 



helpl = 0; 
help2 = 0; 
noconn 1 = 1; 
noconn 2 = 1; 

for ( row=l ; row<glrow; ++row) { 
j = (glcol-1 )*10 + row; 
if ( keygridC j ]==1 ) { 

helpl = helpl + nr_rowsCrow] + noconnl; 
noconnl = 0; 



} 

for 



} 

for 



else noconnl = 1; 

if ( keygridC j ]==1 11 keygridC j+10]==l ) { 

help2 = help2 + 2*nr_rowsC row] + noconn2; 
noconn2 = 0; 

} 

else noconn2 = 1; 



( row=glrow; row<=lastrow; ++row) { 

j = (glcol-1) *10 + row; 

if ( keygridC j ]==1 11 keygridC j-10]==l ) { 

helpl = helpl + 2*nr_rowsC row] + noconnl; 
noconnl = 0 ; 

} 

else noconnl = 1; 
if ( keygridC j ]==1 ) I 

help2 = help2 + nr_rowsCrow] + noconn2; 
noconn2 = 0; 

} 

else noconn2 = 1; 



( subcol =1 ; subcol <=nr_col sC glcol ];++ subcol ) { 

firstnode = subcol_f irstnode ; 
subelement = subcol subelement; 
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coljumpl = key_col jurapl ; 

coljump2 = key_col jump2 ; 

if ( subcol == 1 ) coljumpl = helpl; 

if ( subcol == nr_cols[glcol] ) coljump2 = help2; 

for ( subrow=l ; subrow<=nr_rows[glrow] ;++subrow) { 

cur get ( &row , &col ) ; 

curlocat { 5 , 0 ) ; , 

printfC glcol glrow subcol subrow ); 

curlocat{6,0) ; 

printf { "%4u%8u%8u%9u", glcol, glrow, subcol , subrow) ; 
curlocat ( row, col ) ; 

gl nodenr_of_loc_node[l] = firstnode ; 

gl“nodenr of_loc_node[2] = firstnode + 1; 

gl”nodenr of loc_node[3] = firstnode + 2; 

ql"nodenr of loc node[4] = firstnode + coljumpl +1? 

gl"nodenr"of”loc“node[5] = firstnode + coljumpl + col:jump2 + 2; 

ql“nodenr"of"loc”node[6] = firstnode + coljumpl + col]ump2 + 1; 

gl nodenr of loc_node[7] = firstnode + coljumpl + coljump2; 

gl“nodenr”of loc_node[8] = firstnode + coljumpl; 

iTi 0 rno 2 ry 2 “ nodsnir of loc nodcC^J; 

if{ gl nodenr_of_loc_nodeC5] >= maxnodenrs ) 
maxnodenrs = gl_nodenr of loc nodeCSl; 



/* computation of the coordinates */ 

for (k=l ;k<=8 ;++k) { , n 

X CO keyelemCk] = x_co_keyC { nr_of_keyelement-l ) *8+k] ; 
y”co“keyelem[k] = y_co_key C { nr_of_keyelement-l ) *8+k] ; 



ksi[ 3 ] = - 1.0 + ( 2 . 0 *subcol)/nr_cols[glcol]; 
ksiCl] = ksi[3] - 2 . 0 /nr_cols[glcol ] ; 
ksi[2] = 0.5 * (ksi[l]+ksi[3] ) ; 

eta[3] = 1.0 - ( 2 . 0 *subrow)/nr_rowsCglrow] ; 

etaCl] = eta[3] + 2 . 0 /nr_rows[glrow] ; 
eta[2] = 0.5 * ( etaC 1 ]+eta[ 3 ] ) ; 



for (k=l ; k<=3 ; ++k) 



quad 2 d{ksiCl],eta[k],x_co_keyelem,y_co_keyelem,&xCk],&y[k]) ; 

quad_2d(ksi[2] , eta[3] , x_co_keyelem, y_co_keyelem, Scx[4] , &yL4 J) ; 



for {k=3;k>=l; — k) 



quad ■ 2d(ksiC 3] , eta[k] , x_co_keyelem, y_co_keyelera, &xp^-k] , 
quad"2d(ksi[ 2] , eta[l ] , x_co_keyelem, y_co_keyelem, Scx[3] , &yLa J ) ; 



locCsubelement] = countl + lO0*count2; 

connecC (countl-l)*12 +1] = subelement; 

connecC (countl-l)*12 +2] = idnumber; 

connecC (countl-l)*12 +3] = dummy; 
connecC (countl-1 )*12 +4] = dummy; 

for ( k=l ; k<=8 ; ++k ) { ^ i 

connecC (countl-1) *12 +4 +k] = gl nodenr_of_loc_node 

coordC ( countl- 1 )* 16+2 *k-l ] - xCkJ; 



Ck] 



I 
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coord [ { count 1-1 ) * 16+2*k ] 

} 



= y[k]; 



if { { countl==100 ) 1 1 ( subelement==toteletns ) ) { 

fwrite( (unsigned short *) connec+1 , sizeof { unsigned short), 

12*countl , connf ile ) ; 

fwrite( ( float * )coord+l , sizeof ( float) , 16*countl , coord file) ; 

countl = 0; 

count2 = count2 + 1 ; 

} 

countl = countl + 1; 



coljumpl = coljumpl - 1 

coljump2 = coljump2 + 1 

firstnode = firstnode + 2 

subelement = subelement + 1 

} /* end subrow */ 



subcol_f irstnode = subcol_f ir stnode + coljumpl + coljump2; 
subcol_subelement = subcol_subelement + elemjump; 

} /* end subcol */ 



key_subelement = 
key_coljumpl = 
key coljump2 = 



key_subelement + 
key_coljumpl 
key coljump2 + 



nr_rows[glrow] ; 
nr_rows[ glrow] ; 
nr_rows[glrow] ; 



key_f irstnode = key_f irstnode + 2*nr_rows[glrow] ; 
noconnl = 0; 

memoryl = 0 ; 

for ( row=glrow+l ; row<=lastrow; ++row) { 

if (keygridC { glcol-1 ) *10+row] == 1 ) break; 
if ( glcol == 1 ) { 

memoryl = 1; 

key_f irstnode = key_f irstnode + 1; 

break ; 

} 

if { keygridC ( glcol-2 )* 10+row] == 1 ) { 

memoryl = 1 ; 

key_f irstnode = key_f irstnode + 2*nr_rows [ row] + noconnl; 
noconnl = 0; 

1 



else noconnl = 1; 

} /* end row */ 

if (memoryl == 0) key firstnode = key_f irstnode + noconnl; 



} /* end if there is a keyelement in this glrow in this glcol */ 
} /* end keyelement */ 

glcol_subelement = glcol_subelement + nr_of_subelements[glcol ] ; 

helpl = 0; 

noconnect = 1 ; 

memoryl = 1; 

row = 1 ; 

while ( keygridC glcol *10 + row] == 0 ) { 
if ( keygridC ( glcol-1 )* 10 + row] == 1 ) { 
memoryl = 0; 
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helpl = helpl + 2*nr_rows[row] + noconnect; 

noconnect = 0; 

} /* end if */ 

else noconnect = 1; 
row = row + 1 ; 

} /* end while */ 

if {memoryl==0 ) glcol_f irstnode = meraory2 + helpl + noconnect; 
if (memoryl==l ) glcol_f irstnode = inemory2 + 1; 

} /* end globalcolumn */ 



f seek { connf ile , 0L, 0 ) ; 

fwrite {( unsigned short * ) &maxnodenrs , sizeof ( uns igned short), 1, 

connf ile ) 



f close ( connfile ) ; 

f close { coordf ile ) ; 

locf ile=fopen( natneloc, "w+b" ) ; 

fwrite {{ unsigned short *) loc+1 , sizeof ( unsigned 
f close ( locfile) ; 



short ) , totelems , 
locfile ) 



\n" , s_conn ) ; 

\n" , s_coord) ; 

\n" , s_loc ) ; 

saved on your actual discdrive I \n\n") 



free ( connec ) ; 
free ( coord ) ; 
free( loc) ; 
free (keygrid ) ; 
free { x_co_key ) ; 
free ( y co key) ; 



curlocat(9,0) ; 
printf("The files %s 

printf(" and %s 

printf{" and %s 

printf("are created and 

curlocat (25,0); 



} /* end meshgen */ 
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(8.6.1) mshiptrd.c 

tinclude <stdio.h> 
tinclude <malloc.'h> 
static char s[l5]; 
static char *name = s; 



meshinputread ( s_loc , s_conn, s_coord , x_co_key , y_co_key , keygrid , 
nr_rows , nr_cols , totelems , nrof subelem, lastrow, lasted ) 
char s_loc[],s conn[ ] , s_coord[ ] ; 

float x_co_key[T/ y_co_key [ ] ; 

unsigned short keygrid [ ] , nr_rows[ ] , nr_cols [ ] , * totelems , nrof subelem[ ] 

, * lastrow, *lastcol; 



{ 



FILE * input; 

unsigned short i , j , k , lastrowl , row, col , nr_keyelement s , nr , 

idnumber , maxidnumber=l ; 



int *mess5; 



mess5 =(int * ) calloc ( 540 , sizeof ( int ) ) ; 
if ( mess5==NULL ) 

{ 

printf ( "*calloc ( ) fails for mess5 in mshiptrd.c !"); 
return ; 

} 



wndf rdsk ( "window . 5 " , mess 5 ) ; 
rstwindo (5, 18, 16,62,mess5 ) ; 
curlocat (16,57); 
do { 

scanf ( "%s" , s ) ; 
i=0; 

while (s[i] 1= ’ \0 ' ) ( 

if (s[i] 1= ) { i=i+l ; continue ; } 

break; 

} 

if ( (s[i + l]l = 'k’ ) 11 (s[i + 2]! = ‘e‘ ) I I (s[i+3]l = 'y* ) 1 1 

(filexist(name)==0) ) 

{ 

pr intf ( " \nThe file %s does not exist or has v;rong "); 
printf(" file extens ion 1 \n" , s ) ; 
printf ("Try again : "); 

} 

else break; 

) wh i 1 e ( 1 ) ; 

for ( j=0; j<=i;++j) { 



s loc[j] 


= sCj] 


s conn[j] 


= s[ j] 


s coordtj] 
)" 


= s[ j] 


s loc[i+l] 


= ' 1 ■ 


s loc[i+2] 


= ‘o’ 


s loc[i+3] 


= ‘c’ 


s loc[i+4] 


= ■ \0‘ 
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s connCi+1] 


= 


•c* 


s connCi+2] 


= 


•o' 


s connCi+3] 


= 


' n ' 


s connCi+4] 


_ 1 


\0 ' 


s_coordCi+l ] 


= 


•c* 


s_coordC i+2] 


= 


•o’ 


s coordCi+3] 


= 


•o’ 


s coordCi+4] 


_ 1 


\0’ 


read inputdataf ile 



input = f open ( name, " r” ) ? 
fscanf ( input, "%u" , &nr_keyelement s ) ; 
if (8*nr keyelements >= 399) 

{ 

printf ( "space for x_co_key and y_co_key too small I"); 
return; 

} 

for ( nr=l ; nr <=nr_keyelements ; ++nr ) { 

fscanf ( input, "%u" , Scidnumber ) ; 
keygridC idnumber ] = 1; 

if ( idnumber >= maxidnuraber ) maxidnumber = idnumber; 

} 

fscanf ( input , " %u " , las trow) ; 
fscanf ( input, "%u" , lasted ) ; 
lastrowl = maxidnumber - ( *lastcol-l ) *10; 

/* necessary for later displyay*/ 
for ( row=l ; rov/<= ( * las trow) ; ++row) f scanf ( input , "%u" , &nr_rowsC row] ) 
for (col=l;col<=(*lastcol) ; ++col ) fscanf ( input , " %u" , Scnr_col s[ col] ) 
for { nr=l ; nr <=nr_keyelements ; ++nr ) { • 

for (k=l ;k<=8 ; ++k ) f scanf ( input , "% f " , &x_co_key[ ( nr-1 ) *8+k] ) ; 
for (k=l ;k<=8;++k) f scanf { input ,"% f" , 6cy_co_key[ ( nr-1 ) *3+k] ) ; 
) /* end nr */ 
f close ( input ) ; 

/* compute number of subelements[ ] */ 
for ( col=l ; col <= ( *lastcol ) ; ++col ) { 

nrof subelemC col ] = 0; 

for ( row=l ; row<=( *lastrow) ; ++row) { 
j = (col-l)*10 +row; 
if ( keygridC j] ==1 ) 

nrof subelem[col ] = nrof subelem[col ] + nr_rows[row]*nr_cols[col] ; 

} /* end row */ 

} /* end col */ 

*totelems =0; i 

/* compute toal number of (sub-) elements */ 
for ( col=l ; col <= ( *lastcol ) ;++col) 

*totelems = *totelems + nrof subelemC col ] ; 
clrscrn ( ) ; 

printf ("The computations will be finished when \n\n“); 

printf (" glcol glrow subcol subrow \n"); 

printf( "%4u%8u%8u%9u", *lastcol, lastrowl , nr_colsC *lastcol], ^ 

nr_rowsC lastrowl ]) ; 

curlocat ( 25 , 0 ) ; | 

free (mess 5 ) ; I 
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(8.6.2) quad2d.c 

quad_2d ( ksi , eta , x , y , x0 , y0 ) 
float ksi , eta, x[ ] , y[ ] , *x0 , *y0 ; 
{ 

float n[9]; 

unsigned short function nr; 



*X0 

*y0 

nCl] 


= 


0; 

0; 

0.25 


★ 


( 1-ksi ) 


★ 


( 1+eta ) * 


( eta-ksi-1 ) 


nC2] 


= 


0.50 


★ 


( l-eta*eta) 


★ 


( 1-ksi ) ; 




n[3] 


= 


0.25 


★ 


( 1-ksi ) 


★ 


( 1 -eta ) * 


( -eta-ks i-1 ) 


n[4] 


= 


0.50 


★ 


( l-ksi*ksi ) 


★ 


( 1-eta) ; 




n[5] 


= 


0.25 


★ 


( 1+ksi ) 


★ 


(1-eta) * 


( -eta+ksi-1 ) 


n[6] 


= 


0.50 


★ 


{ l-eta*eta) 


★ 


( 1+ksi ) ; 




n[7] 


= 


0.25 


★ 


( 1+ksi ) 


★ 


(1+eta) * 


( eta+ksi-1 ) 


nC8] 


= 


0.50 


★ 


( l-ksi*ksi ) 


★ 


( 1+eta ) ; 




for 


( function 


nr=l ; function 


nr<=8;++function nr) { 



*x 0 = *x 0 + x[ f unction_nr ] *n[ function_nr7; 
*y0 = *y0 + yCfunction nr ] *n[ function nr]; 
} 

} /* end quad 2d */ 
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(8.7) meshirage.c 

tinclude <stdio.h> 

(rinclude <malloc.h> 

static char s[ 15] , s_loc[ 15 ] , s_conn[ 15 ] , s_coord[ 15 ] ; 
static char *name = s, *nameloc=s__loc, *nameconn = s_conn, 
*naraecoord = s_coord; 
static unsigned short *loc, *connec ; 
static float *coord; 

meshimge { ) 

{ 

FILE *locfile, *connfile, *coordfile; 

unsigned j ,k, countl=0, count2=0, nrsubelein, nodenrs, totelems; 
int i, exitf lag=l ; 
long help; 

connec = (unsigned short *) calloc ( 1201 , sizeof ( unsigned short)); 
coord = (float *) calloc ( 1601 , sizeof ( float )) ; 



clrscrn ( ) ; 

printf("This subroutine displays connectivity 

printf(" ■ and coordinates 

printf ( "after running the meshgenerator . 
printf ( "Enter the name of the meshgenerator input 
printf ("with file extension '.key' : "); 
labell : scanf ( "%s" , s ) ; 



\n") ; 
\n"); 
\n\n") ; 
f ilename\n " ) 



i=0 ; 

wh ile (s[i] != '\0') { 

if (s[i] 1= '.' ) { i=i+l ; continue ; } 

break; 

} 

if((s[i+l] i= 'k' ) 1 1 (s[i+2]!='e' ) 1 1 (s[i+3]l-'y' ) ) { 

printf ("The inputfile has a non valid file extension\n" ) ; 
printf ("Try again : “); 

goto labell; 

} 



for ( j=0; j<=i;++ j ) I 
s_conn[j] = s[j] 
s_coord[j] = s[j] 
s loc[j] = s[j] 
}" 



s conn[i+l] 


= 


* C 


s conn[i+2] 


=: 


* O 


s conn[i+3] 


= 


* n 


s conn[i+4] 


( 


'\o 


s_coord[ i+1 ] 


= 


* c 


s coord[i+2] 


= 


* o 


s coord[i+3] 


= 


* o 


s coord[i+4] 


_ ( 


•\o 


s loc[i+l] 


= 


* 1 


s loc[i+2] 


= 


* o 
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s_l0c[ i+3 ] = ’ c ' ; 
s locC i+4] = ' \0 ' ; 



if { filexist ( nameloc )==0 II filex 

filex 

printf("\n At least on of the 
printf ( " 
printf { " 

printf ("does not exist 1 
printf ("Try again : "); 

goto labell; 

} 



St ( nameconn ) ==0 || 

St ( namecoord )==0 ) { 



files 


%s 


\n" , s conn ) ; 


and 


%s 


\n" , s coord ) ; 


and 


%s 


\n",s loc); 
\n") ; 



locfile = f open ( nameloc r+b" ) ; 
connfile = f open ( nameconn r+b" ) ; 
coordfile = fopen ( namecoord , "r+b" ) ; 

fread (( unsigned short *) &nodenrs , sizeof ( unsigned short), 1, 

connfile ) ; 

fread {( unsigned short *) Sctotelems , sizeof { unsigned short ), 1 , 

connfile ) ; 

loc =(unsigned short * )calloc( totelems+1 , sizeof ( unsigned short)) 
if ( loc==NULL I I connec==NULL I I coord==NULL ) { 

print f {" *calloc ( ) fails in meshimge 1"); 
return; 

) 



fread (( unsigned short *) loc+1 , sizeof ( unsigned short ), totelems , 

locfile ) ; 

for (k=l ;k< = totel :ms ; ++k ) print f (" loc[ %u] : %u \n" , k , loc[k ] ) ; 

getkey { &k ) ; 
clrscrn ( ) ; 

printf(" connectivity image of %s \n\n",s); 

printf ("nr of I belongs to I global nodenumber"); 

printf {" of local node \n"K* 

printf {" subeleraemt I keyelement I "); 

printf (" 1234567 8\n"); 

printf ( "\n" ) ; 

whi le ( ( totelems-count 1 * 1 00 ) >= 100) { 

fread {{ unsigned short *) connec+1 , sizeof { unsigned short), 

1200, connfile ) ; 



for ( j=l ; j <=100 ; ++ j ) ( 

k = (j-l)*12; 

printf{"%6u%12u%12u%5u%5u%5u%5uo5u%5u%5u\n",connec[k+l], 
connec[k+2] , connec[k+5] , connec[k+6] , connec[k+7 ] , connec[k+3] , 
connec[k+9 ] , connec[k+10 ] , connec[k+l 1 ] , connec[k+l 2 ] ) ; 

} 

countl = countl +1; 

1 /* end while */ 

if ( ( totelems-countl *100 ) != 0) ( 

fread (( unsigned short *) connec+1 , sizeof ( uns igned short), 

12* { totelems- 100* countl ) , connfile ) ; 

for ( j=l ; j <= ( totelems-100*countl ) ; ++ j ) { 
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k = (j-l)*12; 

printf ( "%6u%12u%12u%5u%5u%5u%5u%5u%5u%5u\n" , connec[k+l ] , 
connec[k+2] , connec[k+5] , connec[k+6 ] , connec[k+7] , connec[k+8] , 
connecCk+9] , connec[k+10] , connec[k+ll] , connec[k+12] ) ; 

} 

} /* end if */ 

printf ("Hit any key to continue ")? 

getkey (Sc j ) ; 

/* output of coordinates */ 
do { 

clrscrn ( ) ; 

printf ("If you want to see the coordinates of 
printf (" all subelements 

printf (" a specific subelement 

printf (" Exit to mesh mainmenu 
printf("\n Hit (1-3) : "); 

getint ( l,l,-l,l,2,Sci,0,l,2); 
clrscrn ( ) ; 

switch (i) { 

case 1 ; countl = 0; 

while (( totelems-countl * 100 ) >= 100) [ 

f read (( float * ) coord+1 , sizeof ( float ) , 1600 , coordf i le ) ; 

for ( j=l ; j <=100 ; ++ j ) { 

help = ( (countl*100+( j-l)*12)+2L)*2L; 

fseek( connf ile,help, 0 ) r 

f read (( unsigned short *) Scnrsubelem, sizeof ( unsigned short), 1, 

connfile) ; 

printf (" subelem : %4u \n" , nrsubelem) ; 
for ( k=l ; k<=8 ; ++k ) 

printf (" x%lu : %10f y%lu : %10f \n",k, 

coord C ( j-1 ) *16+2*k-l ] ,k , coord [ ( j-1 ) *16+2*k] ) ; 
printf ( "\n" ) ; 

} /* end j */ 

countl = countl +1 ; 

} /* end while */ 

i f (( totelems-countl *100 ) 1= 0) { 

f read ( ( float * ) coord+1 , sizeof ( float ), 16* ( totelems-countl *100 ) , 

coordf ile ) ; 

for ( j=l ; j <= ( totelems -countl*100 ) ; ++ j ) { 

help =( (countl*100+( j-l)*12)+2L)*2L; 

f seek ( connfile , help, 0 ) ; 

f read (( unsigned short *) Scnrsubelem, sizeof ( unsigned short), 1, 

connfile ) ; | 

print f (" subelem : %4u \n" , nrsubelem ) ; 

for (k=l ;k<=8;++k) j 

printf (" x%lu : %10f y?,lu : '510f \n",k, j 

coord [(j-l)*16+2*k-l],k, coord [ (j-l)*16+2*k]); 

printf ( "\n" ) ; , 



HIT\n" ) ; 

1 \n"); 

2 \n\n") 

3 \n"); 
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} /* end j */ 

} /* end if */ 
curlocat (24,0); 

printf("Hit any key to continue 
getkeyU j ) ; 
break ; 

case 2 : do { 

print f { "Enter the number of subelement : "); 
clearkbd ( ) ; 

getint (4 , l,-l,l,2,&i,0,l,9999) ; 

if (i>totelems) printf("\n inputnumber is too large 
printf("\n"); 

} while( i>totelems) ; 



printf ( "\n" ) ; 

f seek ( coordf ile , (loc[i]-l)*64L,0) ; 

f read { { float * ) coord+1 , sizeof ( float ), 16 , coordf ile ) ; 
for (k=l ;k<=8;++k) 

printf(" x%lu : %10f y%lu ; %10f \n",k 

coord [ 2*k-l ] ,k, coord [ 2*k] ) 



printf ("\n"); 



curlocat (24,0) ; 

printf ("Hit any key to continue .."); 
getkey( & j ) ; 
break ; 

case 3 : exitflag = 0; 

default : break; 



} /* end switch */ 



) while (exitflag) ; 



fclose ( locf i le ) ; 
fclose(coordf ile) ; 
fclose ( connf ile ) ; 
free ( connec ) ; 
free ( coord ) ; 
free ( loc ) ; 

) 
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(8.8) meshgraf.c 

#include <stdio.h> 

#include <malloc.h> 

static char s[ 15 ] , s_conn[ 15 ] , s_coord[ 1 5 ] ; 

static char *name = s , *nameconn=s_conn , *namecoord=s_coord ; 
static float *coord; 

meshgraf ( ) 

{ 

FILE *connf ile , *coordfile; 

unsigned short totelems , max_gl_nodenr ; 

unsigned short countl=0 , k , 1 , j ; 

int i,x[8],yC8], factor , f actor_x , factor_y r 

float min_x=1000. ,min_y=1000 . ,max_x=-1000 . , max_y=-10O0 . ,help; 

coord = (float *) calloc ( 1601 , sizeof ( float )) ; 
if ( coord==NULL ) { 

printf ( "*calloc ( ) fails in meshgraf I"); 
return? 

} 



clrscrn ( ) ; 

printf ("This subroutine does create a graphic display"); 
printf (" of the mesh \n"); 

printf ( "generated by the meshgenerator ! \n\n"); 

printf ( "Enter the name of the meshgenerator input filename\n" ) 
printf (" ( file extension '.key' ) ; "); 

labell : scanf ( "%s " , s ) ; 

i=0 ; 

while (s[i] 1= ' \0 ' ) { 

if (s[i] != '.' ) { i=i+l ? continue ; } 

break; 

} 

if((sCi + l] 1= 'k' ) I I (s[i+2]l = 'e' ) I I (s[i+3]l = 'y' ) ) { 

printf ("The inputfile has a non valid file extension\n" ) ; 
printf ("Try again ; "); 
goto labell; 

} 

for ( j=0; j<=i;++ j ) { 

s_conn[j] = sCj]; 
s coord[j] = s[j]; 

}" 

s_conn[i+l] = 'c'; 
s_conn[i+2] = 'o'; 
s_conn[i+3] = 'n'; 
s_conn[i+4] ='\0'; 
s_coord[i+l] = 'c'; 
s_coord[i+2] = 'o'; 
s_coord[i+3] = 'o'; 
s_coord[i+4] ='\0'; 

if ( f i lexist ( nameconn ) == 0 || f ilexist ( namecoord ) == 0) { 
printf (" \nEither %s \n",s_conn); 
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\n" , s coord ) ; 



\n") ; 



printf ( "or %s 

printf("does not exist I 
printf ("Try again : "); 

goto labell; 

} 



connfile = fopen ( nameconn , " r+b" ) ; 
coordfile = fopen(namecoord,"r+b"); 
f seek ( connfile , 2L, 0 ) ; 

fread (( unsigned short *) &totelems , sizeof ( unsigned short),l, 

connfile 



k = 8*totelems +1; 

while (( totelems - countl*100) >= 100 ) { 

fread (( float * ) coord+l , sizeof ( float) , 1600, coordf i le ) ; 
for ( k=l ; k<=800 ? ++k ) { 



if ( min_x > = 
if ( min_y >= 
if ( max_x < = 
if ( max y <= 

) " 

countl = countl 



coord[ 2*k-l ] 
coord[2*k] 
coord [ 2*k-l ] 
coordC 2*k] 

+ 1 ; 



) min_x = coordC 2 *k-l ] ; 
) min_y = coord[2*k]; 

) max_x = coordC 2 *k-l ] : 
) max y = coord[2*k]; 



if ((totelems - 100*countl) 1=0 ) { 

f read ( ( float * ) coord+l ,sizeof(float) ,16*( totelems- 100* countl ) , 

coordfile) ; 

for (k=l ; k<=8* ( totelems- 100* countl ) 



if 

if 

if 

if 



min 
min 
max 
max 
} ■ 



X >= coord[2*k-l] 
y >= coord[2*k] 

X <= coord[2*k-l] 
y <= coord[2*k] 



++k) { 

) min_x = coordC 2*k-l ] ; 
) min_y = coordC2*k]; 

) max_x = coordC 2 *k-l ] ; 
) max y = coordC 2*k]; 



factor_x = 303 . 0/ ( rnax_x-min_x ) ; 

factor_y = 183 .0/ (max_y-min_y ) ; 

if( factor_x >= factor_y) factor = factor_y; 

else factor = factor_x; 

countl = 0; 

grafinit(2,0,0) ; 

f seek ( coordfile, 0L , 0 ) ; 

while (( totelems-100*countl ) >= 109 ){ 

fread ( ( float * ) coord+l , sizeof ( float ) , 1600 , coord file) ; 
for (k=l ; k<=lO0 ; ++k ) { 

for (1 = 1;K=8;++1) { 

xCl-l] = 8 + factor *( coordC (k-l)*16 + 2*1-1 ] - min x ) ; 
yCl-1] = 191 - factor * (coordC (k-l)*16+2*l ] - min yT ; 

} 

grpoly(x,y,8, 1 ) ; 

} 

countl = countl + 1; 

} /* end while */ 

if ( ( totelems-100*countl ) != O ){ 

fread (( float * ) coord+l , sizeof ( float ), 16* ( tote lems-1 GO* count 1 ) 

coordfile ) 
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for (k=l ; k<=totelems-100*countl ; ++k ) { 

for (1=1;K=8;++1) { 

x[l-l] = 8 + factor *( coord[ (k-l)*16 + 2*1-1 
yCl"l] = 191 - factor * (coord[ (k-l)*16+2*l ] 

} 

grpoly(x,y,8,l) ; 

} 

} /* end if*/ 

grxlab("Hit any key to continue ..",7,54,0); 
i = getkey( & j ) ; 
graf init (0,0,0); 
fclose ( connf ile ) ; 
f close ( coordfile ) ; 

free ( coord ) ; 

} 



] - min X ) 
- min yT; 
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